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Abstract 

Background: Studies of oyster microbiomes have revealed that a limited number of microbes, including 
pathogens, can dominate microbial communities in host tissues such as gills and gut. Much of the bacterial 
diversity however remains underexplored and unexplained, although environmental conditions and host genetics 
have been implicated. We used 454 next generation 16S rRNA amplicon sequencing of individually tagged PCR 
reactions to explore the diversity of bacterial communities in gill tissue of the invasive Pacific oyster Crassostrea 
gigas stemming from genetically differentiated beds under ambient outdoor conditions and after a multifaceted 
disturbance treatment imposing stress on the host. 

Results: While the gill associated microbial communities in oysters were dominated by few abundant taxa 
(i.e. Sphingomonas, Mycoplasma) the distribution of rare bacterial groups correlated to relatedness between the 
hosts under ambient conditions. Exposing the host to disturbance broke apart this relationship by removing 
rare phylotypes thereby reducing overall microbial diversity. Shifts in the microbiome composition in response 
to stress did not result in a net increase in genera known to contain potentially pathogenic strains. 

Conclusion: The decrease in microbial diversity and the disassociation between population genetic structure of 
the hosts and their associated microbiome suggest that disturbance (i.e. stress) may play a significant role for 
the assembly of the natural microbiome. Such community shifts may in turn also feed back on the course of 
disease and the occurrence of mass mortality events in oyster populations. 

Keywords: Microbiota, Population structure, Stress, Pathogen, Biological invasion, Pacific oyster, Crassostrea 
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Background 

Emerging diseases of marine organisms often manifest 
in mass mortalities associated with environmental per- 
turbations such as heat stress events [1]. This also ap- 
plies to marine bivalves where infectious agents cause 
detrimental effects by profiting from increased tempera- 
tures in combination with a weakened immune response 
of the host [2,3]. Prominent examples for such mass 
mortalities are summer mortalities' of farmed and wild 
Pacific oysters Crassostrea gigas in several localities 
worldwide [4-6]. Here, the outcome of an infection is 
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thought to be driven by a complex interplay of abiotic 
factors (e.g., temperature) and biotic factors (e.g., host 
genetic or immune system effects [7,8] and/or repro- 
ductive state [9]). More recently, host-associated micro- 
biota have also been suggested to play an important role 
in determining host fitness [10,11]. Such effects can be 
mediated by providing additional energy sources by 
chemosynthesis [12] but also in defence against disease 
by either preventing establishment of pathogens or dir- 
ectly attacking them with antimicrobial effector mole- 
cules [13]. The use of probiotics in bivalve aquaculture 
has therefore been discussed as a means of preventing 
loss due to disease [14]. However, relatively little is 
known about microbial communities of native popula- 
tions and their response to environmental perturbations. 
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Microbial communities residing in different organs of 
several oyster species have only recently been described 
by using molecular, culture independent techniques 
[15-17] that allow intra- and interspecies comparisons 
[18] and the exploration of environmental factors, such 
as temperature [19]. For example, oysters invading the 
Mediterranean from the Indian ocean maintained some 
of their associated microbes throughout the invasion 
process [18]. This is not self evident as biological inva- 
sions inevitably represent a drastic shift in abiotic and 
biotic conditions including the exposure to different mi- 
crobes in the new environment. 

Using next generation amplicon sequencing of indi- 
vidually tagged 16S rRNA-PCR reactions [20], we here 
assessed the combined effects of host population and 
disturbance/host stress on the microbial communities 
associated with gill tissue of Pacific oysters Crassostrea 
gigas stemming from populations only very recently in- 
vading the North Sea. The invasion of Pacific oysters 
into the Wadden Sea part of the North Sea originated 
from aquaculture activities in the 1990s [21], and today 
Pacific oysters locally represent the dominant epibenthic 
bivalve species [22]. Oyster populations in the northern 
and southern parts of the Wadden Sea stem from two 
genetically distinct invasion sources [23]. These separate 
invasions are also interesting in terms of summer mor- 
tality events because summer mortality has been ob- 
served only in southern populations so far [24] . 

Individual microbial communities can also be influenced 
by host genetics, either between populations (i.e. phylo- 
geography and genetic differentiation) [25] or within pop- 
ulations (i.e. relatedness). Strong skew in reproductive 
success among individual breeders [26] is common in 
marine bivalves displaying high juvenile mortality (i.e. type 
III survivor curves) and can lead to increased genetic dif- 
ferentiation. In turn this can also lead to genetic differenti- 
ation on small spatial scales and therefore we here 
compare microbial communities in oysters from differ- 
ent reefs that most likely originated from different 
spatfall events. Our sampling scheme allowed us to 
evaluate the relative importance of host population 
genetic structure independent of confounding effects 
of geography. We investigated a total of 40 individual 
oyster microbiomes within three separate oyster reefs 
stemming from two tidal basins in the northern Wadden 
Sea. By exposing half of the oysters to a disturbance treat- 
ment, we tested if stress in combination with environmen- 
tal change causes a shift in the microbial communities 
and if such a shift is associated with an increase in the 
abundances of potentially pathogenic bacteria during pe- 
riods of stress. This could potentially reveal whether mor- 
tality events originate from environmental or intrinsic 
reservoirs and if such events are possibly associated with 
the demise of beneficial microbes. 



The artificially induced microbial community shift can 
thus be used to compare reaction norms of microbial 
communities in naturally replicated host genotypes 
across genetically differentiated host populations. Our 
detailed objectives were 1) to test the differentiation of 
individual host-associated microbial communities ac- 
cording to population and individual genetic differenti- 
ation (i.e.phylogeography and relatedness, respectively) 
and 2) to identify bacterial taxa from those communi- 
ties that respond similarly to disturbance stress across 
individuals and host populations in order to evaluate 
whether such pulsed perturbations suffice to promote a 
rise in opportunistic and potentially pathogenic strains. 

Methods 

Sampling & experimental procedures 

To explore the relationship between host genetic differ- 
entiation and microbiome composition in response to 
environmental stress we collected oysters on 18 th and 
23 rd of January 2008 from three oyster beds in the 
northern Wadden Sea covering two tidal basins, the 
Sylt-Romo-Bight (Diedrichsenbank - DB 55° 02' 32.13" 
N, 08° 27' 02.86" E, Oddewatt OW 55° 01' 41.20" N, 
08° 26' 17.31" E) and the Hornum Deep (Puan Klent 
PK 54° 47' 29.59" N, 08° 18' 18.52" E, see Figure 1). 
We chose to collect oysters in winter because diversity 
and abundance of pathogenic strains are correlated with 
temperature [27] and the input of transient open water 
pathogens could potentially be minimised this way. 
From each bed we collected 20 oysters by picking single, 
unattached individuals from soft-bottom mud flats. After 
collection half of the oysters were frozen (-20°C) while 
the other half was transferred to large buckets (20 L) 
filled with sand-filtered seawater (salinity 29%o). We kept 
groups of oysters in these buckets under constant aer- 
ation at densities of 10 oyster /bucket. To minimise al- 
lochthonous input of microbes and facilitate spread of 
potential pathogens we decided to use static conditions 
with no flow-through and did not feed the oysters during 
the experimental treatment. All experimental animals 
were exposed to a heat-shock treatment by increasing 
water temperature from ambient 2°C to 26°C over a time 
span of 10 days, before individuals were frozen at -20°C. 
We chose this steep temperature increase to maximise 
heat-induced stress for the host and to allow potential 
pathogens to proliferate since temperatures of >20° are 
often associated with pathogen induced mass mortalities 
[24,28]. Our disturbance treatment thus combined aspects 
of transfer, food and heat stress. All experiments complied 
with German legal standards. For genetic analyses a small 
piece of gill tissue was removed from each individual oys- 
ter and DNA was extracted using the Wizard Genomic 
DNA Purification kit (Promega, Mannheim) following the 
manufacturers instructions. We decided to use gill tissue 
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Figure 1 Geographic location and genetic differentiation 
between investigated oyster beds. Stars indicate the location 
of the oyster beds and boxes the pairwise genetic differentiation 
(F S j) between host populations. ** p < 0.001, * p < 0.01, ns: 
not significant. 
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because gills constitute large contact surfaces to the sur- 
rounding water and should thus capture both, resident 
bacteria as well as bacteria from the environment. Fur- 
thermore, it has been shown that gill microbiota of Medi- 
terranean oysters are more distinct from surrounding 
waters than those associated with gut tissue [18]. We used 
14 oysters per bed (7 ambient ones frozen immediately 
and 7 exposed to disturbance treatment in the lab) for 
genetic analysis and microbiome sequencing. 

Genotyping of host microsatellite markers 

To measure host population structure we amplified ten 
polymorphic microsatellite loci [29] covering nine link- 
age groups. These included Cgl57 (LG I 141.5 cM), 
Cg_194 (LG I 28.7 cM), Cgl36 (LG II), Cgl93 (LG III), 
Cg_164 (LG IV), Cgl73 (LG V), Cg_i28 (LGVI), Cgi24 
(LG VII), Cgl72 (LG IX) and Cg_181 (LG X) [30]. Loci 
were amplified using M13-tail labelling [31] in 20 ul 
volume using 4 ul of 5x concentrated buffer containing 
MgCl 2 (Promega, Mannheim), 10 nM of each dNTP, 
13.9 ul of H 2 0, 5 nM of each locus specific primer, 8 pmol 
of one M13-tail labelled with either FAM, VIC, NED or 
PET fluorophores and 1 unit of GoTaq Polymerase 



(Promega, Mannheim). Cycling used a standard protocol 
consisting of 5 min at 94°C followed by 28 cycles at 94°C 
(30 s) / 56°C (45 s) / 72°C (45 s). M13 tail incorporation 
was achieved in 8 cycles at 94°C (30 s) / 53°C (45 s) / 72°C 
(45 s), and a final extension at 72°C for 10 min. PCR prod- 
ucts were pooled into sets of four loci with differently 
coloured labels and separated on a ABI Prism 3130 XL 
(LifeTechnologies, Darmstadt) capillary sequencer using a 
LIZ 500 size standard. Product sizes were scored using the 
GeneMarker vl.90 software (SoftGenetics) and pairwise 
genetic differentiation between populations was calculated 
as Wrights fixation indices (F ST ) according to Weir & 
Cockerham [32]. Pacific oysters are known to harbour 
substantial amounts of null alleles [33] that could bias any 
estimate of population differentiation. We therefore esti- 
mated the frequency of null alleles within our sample 
using the software freeNA [34]. Frequencies of null alleles 
were small for all loci and populations (range: 0 - 0.06) 
except for locus Cgi-194 where estimates were higher 
within all oyster beds (range: 0.05 - 0.15). We therefore 
excluded this locus from the analysis, and only report the 
corrected F S t values after removal of loci with high fre- 
quencies of null alleles in all populations. Genetic distance 
between individuals was calculated as geometric AMOVA 
distances: 

5=1 

where distances between individual genotypes j,k are 
summed over S loci. Calculations were performed using 
the R package Gstudio. 

Amplicon sequencing of microbial communities 

Microbial diversity and composition was measured 
within a standardised amount of genomic DNA (30 ng). 
We used an informative part of the 16S rRNA gene 
spanning position 535-904/912 in the E. coli genome 
covering the variable regions V3 and V4 for ribotyping. 
Using a PCR product of this length increases the preci- 
sion of taxonomic assignment [36] and should provide 
high resolution due to the inclusion of two variable re- 
gions. Initial testing of these primers revealed that they 
preferentially amplified host 18S rRNA (24 out of 24 
randomly picked clones). To avoid cross-amplification 
we performed a nested PCR starting with primers spe- 
cific for Eubacteria (27f: AGA GTT TGA TCA TGG 
CTC AG and 1492r: GGT TAC CTT GTT ACG ACT 
T). The first nested PCR consisted of 30 ng of genomic 
DNA, 0.05 ul of Hot start taq (5 unit/ul, Promega), 
1 mM of each dNTP, 4 ul of reaction buffer (Promega), 
1 ul of each forward and reverse primers (5 uM) and 
11.5 ul of molecular grade water. Cycling started with an 
initial denaturation and hot start activation of 10 min at 
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95°C followed by a low number of 16 cycles of 30 s de- 
naturation at 95°C, 30 s at 50°C and 90 s at 72°C and a 
final extension of 10 min at 72°C. One \A of each PCR 
product was then diluted in 99 \A of molecular grade 
water before the internal stretch was amplified for 
454 sequencing. Here, each individual microbiome 
was tagged by a unique combination of multiplex 
identifiers (MID, Roche, Basel, CH) integrated into 
forward and reverse primers [37,38]. We used a total 
of 20 tagged primers consisting of the Titanium B se- 
quencing adaptor (Roche, Basel), the 454 sequencing 
key, a MID tag and the gene-specific sequence. Hence, 
an example of a forward primer would have the follow- 
ing sequence: 5 -CCATCTCATCCCTGCGTGTCTCC- 
GAC TCAG ACGAGTGCGT CCACGAGCCGCGGT 
AAT -3' and a reverse primer: 5 -CCTATCCCCTGT 
GTGCCTTGGCAGTCTCAG TCAG ACGAGTGCGT 
CCGTCAATTCMTTTAAGTTT-3 ', with the 454 se- 
quencing key in italics, the MID tag in bold and gene 
specific sequence underlined. Combinations of forward 
and reverse MIDs were random with respect to treat- 
ment and oyster bed. Therefore any amplification bias 
introduced by the MID will be randomly distributed 
among groups. 

After amplification single PCR reactions were purified 
using the MinElute 96 kit (Qiagen, Hilden) before 2 (il of 
each elution was used for pooling. To eliminate remaining 
primer-dimer both pools were purified again using Wizard 
PCR clean-up system (Promega, Mannheim) following 
the manufacturer s instructions. After confirming the sole 
presence of the desired PCR product without any traces of 
primer by gel electrophoresis, the pool of individually bar- 
coded PCR reactions were sequenced on the 454 FLX 
genome sequencer (Roche, Basel, CH) using Titanium 
chemistry. Sequencing was performed by GATC Biotech 
(Konstanz, Germany). 

Data analysis 

Assignment of reads to individual PCRs was done using 
modified python scripts from the cogent package. In 
short, within each raw read we looked for the presence of 
both primers ensuring complete sequencing of the PCR 
product. Afterwards, we identified individuals by deter- 
mining combinations of MID tags allowing for a max- 
imum hemming distance of one in each MID tag. After 
correct assignment of single reads to an individual oysters, 
we used the AmpliconNoise pipeline [39] to remove pyro- 
sequencing and PCR noise and Perseus to remove 
chimeric sequences using default parameters except for 
alpha and beta values for false discovery detection in Per- 
seus, which were set to -7.5 and 0.5, respectively. Reads 
were trimmed by cutting off their forward and reverse 
primers. 



We used scripts from the Qiime package [40] for the 
analysis of microbial diversity. In detail, we clustered se- 
quences into operational taxonomic units (OTUs) using 
UCLUST using unique and cleaned reads assigned to in- 
dividual oysters. Taxonomic assignment of OTUs found 
for each individual oyster was done using the naive 
Bayesian Classifier [41]. We used an assignment cer- 
tainty threshold of 60% for each taxonomic classification. 
As singleton reads overestimate the contribution of rare 
phylotypes [42] we removed singleton reads. All analyses 
were then based on the resulting OTU table to account 
for small strain specific differences and was used to cal- 
culate observed bacterial diversity (Shannons H'). Suffi- 
cient sampling of observed diversity was confirmed by 
rarefactions based on group specific microbiomes. Po- 
tentially pathogenic OTUs were identified by genus clas- 
sifications and pooled according to genus affiliation. We 
used previously described genera of pathogenic bacteria 
in shellfish [3] and other marine organisms [43] to iden- 
tify such potentially pathogenic bacteria. These included 
Arcobacter spp., Citrobacter spp., Corynebacterium spp., 
Escherichia spp., Halomonas spp., Micrococcus spp., 
Mycoplasma spp., Photobacterium spp., Pseudoalteromo- 
nas spp., Pseudomonas spp., Shewanella spp., Staphylo- 
coccus spp., Streptococcus spp., Tenacibaculum spp.. We 
used non-metric multidimensional scaling from the 
vegan R package to visualise distance matrices (Horn- 
Morisita distances, Wisconsin double square root trans- 
formation) between individual microbiomes. Statistical 
differences between treatments and oyster beds were 
analysed by means of multivariate permutational ANOVA 
{adonis function, Horn-Morisita distances) and com- 
parisons between distance matrices were based on 
non-parametric Mantel tests or procrustes rotations of 
ordinations. To account for differences in sequencing depth 
between libraries we also resampled all communities to the 
lowest coverage using the perl script daisychopper (avail- 
able at http://www.genomics.ceh.ac.uk/GeneSwytch/Tools. 
html). To further account for differences in library size, 
analyses relying on the abundance of OTUs (e.g. abun- 
dance - occupancy analyses) were based on relative 
abundances of ln-transformed read numbers within each 
oyster. All analyses were performed in R [44]. 

Results 

Host genetic differentiation 

We found significant genetic differentiation (F S t) in two 
out of the three pairwise comparisons between oyster 
beds (Figure 1). Interestingly with a F 5:r value of 0.043 
(P < 0.001) the largest pairwise differentiation was ob- 
served between the two oyster beds found closest to 
each other, i.e. Diedrichsenbank (DB) and Oddewatt 
(OW, geographic distance 2.5 km) while the genetic dif- 
ferentiation to a different tidal basin was lower (OW-PK: 
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¥ S t= 0.026, P = 0.002) or not even significant (DB-PK: 
F 5r = 0.009, P = 0.124, Figure 1). Similar results were ob- 
tained using individual genetic differentiation based on 
AMOVA distances and multivariate Permanova of the 
resulting distance matrix, with genetic variation among 
populations explaining 11.9% of the overall variation 
(P < 0.001 based on 1000 permutations). In concord- 
ance with the expectation of random sampling before 
treatment assignment we found no significant differ- 
ence between "ambient" and "disturbed" oysters in 
terms of their genetic variation (R 2 = 0.031, P = 0.159 
based on 1000 permutations) and no significant inter- 
action effect (R 2 = 0.053, P = 0.257 based on 1000 per- 
mutations). Due to high within locus polymorphism 
the majority of variation was found among individuals 
(R 2 = 0.797). 

Microbial communities of oysters before and 
after disturbance 

Out of the 52,092 reads that could successfully be 
assigned to an amplicon library for each individual, 
38,029 reads passed our quality selection and de-noising 
criteria for further analysis. The resulting average library 
size per individual was 825 ± 80. With a total number of 
4,464 unique operational taxonomic units (OTUs) dis- 
tributed over 213 genera, microbial species richness was 
very high. However, only few OTUs occurred frequently 
and most OTUs occurred rarely (<1% within whole data 
set). After rigorous de-noising of our sequencing data 



we potentially underestimated species richness of the re- 
spective communities, but we could reliably calculate di- 
versity (Shannons H') for most experimental groups 
(Figure 2A). Microbial diversity was significantly lower 
in oysters originating from DB (GLM, F 2 ,36 = 3.55, P = 0.039) 
especially under ambient conditions (Figure 2A,B). The dis- 
turbance treatment led to a significant decrease of bacterial 
diversity in oysters from all beds (Figure 2B, disturbance: 
GLM F 1>36 = 7.52, P = 0.009, disturbance x oyster bed 
interaction: F 2>36 = 0.80, P = 0.456). 

Non-metric multidimensional scaling of the full bac- 
terial communities from individual oysters suggested 
that communities were differentiated by treatment along 
both axes (Figure 3), which was confirmed by Perma- 
nova (effect of disturbance: R 2 = 0.077, P = 0.006). Clus- 
tering of ambient group centroids in the ordination 
suggests that initially there was no significant difference 
between beds and large variation within beds under am- 
bient conditions (Figure 3, Permanova, effect of oyster 
bed: R 2 = 0.058, P = 0.211). Community shifts in response 
to the disturbance treatment however seemed to depend 
on the initial community composition and were differen- 
tiated in the NMDS as oysters from PK shifted vertically, 
while shifts occurred mainly horizontally for oysters 
stemming from OW and DB. This effect was however only 
marginally significant in the overall analysis (Permanova, 
disturbance x oyster bed interaction: R 2 = 0.076, P = 0.073). 
Similar results were obtained with rarefied communities 
(n = 245 reads per library, disturbance effect: R 2 = 0.078, 
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Figure 2 Bacterial diversity (Shannon's H') of oyster gill microbiota stemming from different oyster beds. A) Rarefaction curves of 
Shannon's H' in different oyster beds under ambient field conditions and after disturbance. Shown are rarefied means for treatment and origin 
groups from 10 resamples with a maximum number corresponding to the lowest coverage of a single microbiome in each group. Solid lines 
represent ambient conditions and dashed lines disturbed microbial communities. B) Observed values of Shannon's H' for individual oysters 
stemming from different oyster beds (mean±se) showing significant differences between oyster beds (F 2)36 = 3.55, P = 0.039) and a significant 
decrease of diversity after disturbance (F 136 = 7.52, P = 0.009). 
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Figure 3 Non-metric multidimensional scaling of bacterial 
communities associated with oyster gill tissue. Ordination was 
based on Horn-Morisita distances from unique OTUs after Wisconsin 
double standardisation and square root transformation. Symbols 
show communities of single oysters with circles representing ambient 
communities and triangles representing disturbed communities. Solid 
and dashed lines connect single communities with group centroids. 
Colours code for different oyster beds. 



P = 0.009, oyster bed effect: R 2 = 0.054, P = 0.244, disturb- 
ance x bed interaction: R 2 = 0.076, P = 0.081). 

Proteobacteria represented the numerically most abun- 
dant phylum in both ambient and disturbed conditions 
(Figure 4). Numerical abundance of Proteobacteria was 
owed to the fact that the overwhelming majority of OTUs 
were affiliated with the genus Sphingomonas (30.6 - 64.1% 
for each treatment group, Figure 4) and only few other 
taxa reached comparably high numbers (e.g. Flavobacteria 
{Barter oidetes)). Several taxa were characteristic for spe- 
cific oyster beds or shifts during disturbance treatment 
(Figure 4). Flavobacteria (Barter oidetes), for example, 
were common in OW and PK in ambient conditions rare 
in DB. All beds contained several genera unique to each 
treatment with ambient communities having higher pro- 
portions of unique taxa reflecting higher overall diversity. 
Disturbed communities from DB and OW oysters were 
shaped by OTUs associated with the genus Mycoplasma 
(Tenericutes) while Planctomycetales were characteristic 
for disturbed communities of PK oysters (Figure 4). 

Contrary to our expectation we did not observe a sig- 
nificant increase in the proportion of reads containing 
potentially pathogenic bacterial genera after the disturb- 
ance treatment (paired £-test, t = 0.990, df = 17, P = 
0.336) nor did we find an increase in their taxonomic 
abundance (DB: 2 taxa unique in ambient communities 
vs. 2 taxa in disturbed communities, OW: 4 vs. 2, PK: 7 



vs. 2, Figure 4). While the overall load of genera contain- 
ing known pathogenic strains did not change signifi- 
cantly, single genera increased or decreased strongly in 
response to the disturbance (Figure 4). Reads classified 
as Mycoplasma increased strongly in abundance while 
other well established shellfish pathogens like Vibrio 
were very rare (Figure 4, frequency 0.013%). 

Abundance (i.e., how frequent an OTU occurs in a 
host) is often positively correlated to occupancy (i.e. the 
number of hosts an OTU is observed in) [45]. We found 
such a significant relationship between the mean relative 
abundance of OTUs in single oysters and the number of 
oysters they occurred in (occupancy) only after disturb- 
ance (Spearman's rank correlation: p = 0.175, P < 0.001) 
while ambient bacterial communities did not show such 
a relationship (Spearman's rank correlation: p = -0.004, 
P = 0.931). In both environments we could identify some 
generalist taxa (moderately abundant in more than 50% 
of hosts [46,47]). Specialist taxa (highly abundant in less 
than 25% of hosts) were rare under ambient conditions 
but we could observe a shift towards increased specialisa- 
tion in disturbed communities that was mainly associated 
with a steep increase in relative abundance of OTUs asso- 
ciated to the genus Mycoplasma (Figure 5A). 

Overall, only few OTUs were observed in both treat- 
ments (n = 298 corresponding to 6.7%) and we could ob- 
serve a net increase in relative OTU abundance (paired 
£-test, mean difference = 0.19, t = 3.96, df=297, P< 
0.001) but a net decrease in occupancy (paired £-test, 
mean difference = -0.32, t = -2.19, df = 297, P = 0.029). 
Taxonomic differences between ambient and disturbed 
communities were mainly associated to a decrease of 
Bacteroidetes (especially Flavobacteria) and the in- 
crease of Actinobacteria in disturbed oysters (Figures 4 
and 5B), which corresponded to the disappearance of 
generalist Flavobacteria after disturbance (Figure 5A). 

Correlation of microbial community and host population 
genetic structure 

In contrast to host population structure (Figure 1) we 
did not find a significant difference in microbial commu- 
nity structure on the level of oyster beds (Figure 3). Con- 
sidering that most genetic as well as microbial community 
variation was partitioned between individuals, microbial 
communities could also associate with individual geno- 
types within populations rather than with geographically 
and genetically separated host populations. Accordingly 
we found a significant correlation of individual pairwise 
genetic distances (AMOVA) and microbial community 
distances (Bray-Curtis dissimilarity) for ambient oysters 
using non-parametric Spearman's rank correlation reflect- 
ing the non-normal distribution of microbial community 
distances (Mantel test: R = 0.137, P = 0.045). This result 
was supported by a correlation of symmetric procrustes 
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(See figure on previous page.) 

Figure 4 Association network of bacterial taxa (genus level) in ambient and disturbance treatments of the three different oyster beds 
(DB, OW, PK). Taxa are shown as circles with colour-coded phylogenetic relationship and size reflecting overall relative abundance (ln(x+ 1) 
transformed) from a rarefied resampled data set. Lines indicate the occurrence in the respective treatment. Hence, taxa only related to one 
treatment occurred exclusively in either ambient or disturbed oysters while taxa related to both treatments occurred before and after the disturbance. 
Width of lines corresponds to the proportion of each taxon within each treatment. Red edges indicate significant distribution bias towards one 
treatment group. Hexagons represent genera known to contain pathogenic strains (printed in bold) and numbers below each group give the number 
of genera followed by the number of potential pathogens. 

V ) 



rotations of both ordinations (R = 0.48, P = 0.018 based on 
1000 permutations). Such a result was not observed for 
disturbed oysters (Mantel test: R = -0.07, P = 0.756, Pro- 
crustes rotation R = 0.19, P = 0.714 based on 1000 permu- 
tations) indicating that original communities may have 
adjusted to different host genotypes while these associ- 
ation broke apart as a result of disturbance. 

We subsequently tested whether rare or common 
components of the bacterial communities were respon- 
sible for the observed correlation and removed OTUs 
in a sliding window approach based on their abun- 
dance. In detail, we first removed OTUs that occurred 
only twice in the data set and repeated the correlation 
analysis for both ambient and disturbed oysters. This 
procedure was iterated with increasing abundance cut- 
off values up to an abundance threshold of 100, which 
represents a reasonable upper limit because communi- 
ties contained only few taxa after this procedure and 
only changed little with higher thresholds. We only 



found significant positive correlations for communities 
containing rare OTUs (overall abundance threshold 2-4) 
while all disturbed communities correlated negatively with 
genetic distance among individuals (Figure 6). 

Discussion 

Linkages between microbial communities and 
host genetics 

Using a barcoded 16S rRNA 454 amplicon sequencing 
approach, we showed that the undisturbed ambient mi- 
crobial communities associated with gill tissue of invasive 
Pacific oysters Crassostrea gigas from three genetically dif- 
ferentiated reefs did not follow genetic differentiation of 
the host populations, although environmental variation 
and host population structure have previously been impli- 
cated in the composition of gut microbial communities 
[25]. Oyster gill microbiota, on the other hand, harboured 
a substantial amount of variation between individuals 
(Figures 2 and 3). The between individual variation in 
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Figure 5 Relationships between abundance and occupancy of OTUs recovered from oyster gill tissue. A) Abundance occupancy plot 
showing the relative mean abundance ((In + 1) transformed) of each OTU as a function of occupancy (i.e., from how many oysters it was 
recovered) for ambient (blue circles) and disturbed conditions (red triangles). Filled symbols mark generalists (abundance less than 1% in more 
50% of oysters) and specialist (highly abundant in few oysters) OTUs. Pie charts show the taxonomic affiliation of generalists and specialists, 
where the size of the pie corresponds to the number of OTUs. B) Taxonomic composition of all taxa that increased (upper panel) or decreased 
(lower panel) in abundance and occupancy. Pie size represents number of OTUs found in each group and colours code for different phyla. 



Wegner et al. BMC Microbiology 2013, 13:252 
http://www.biomedcentral.com/1471-2180/13/252 



Page 9 of 12 




0 20 40 60 80 100 



Abundance of OTUs removed 
Figure 6 Correlation coefficients (Spearman's) between genetic 
distance among individuals and similarity of microbial 
communities associated with host gill tissue. The blue and red 
lines represent ambient and disturbed communities, respectively. 
OTUs were iteratively removed with increasing abundance 
thresholds and significance of each correlation was assessed by 
Mantel tests with 1000 randomisations. Significant correlations (p < 
0.05) are shown as triangles and could only be observed for 
correlations containing rare parts of the ambient communities. 



microbial community composition correlated with genetic 
relatedness of the oysters, suggesting that microbial com- 
munities might assemble according to individual hosts or 
even host genotypes. Stable host associations have been 
reported for several gut microbiota in a variety of host 
species [48-51]. The human gut bacterial community, for 
example, is considered to be stable over extended periods 
of time, but is also unique for each individual [51] and 
similar between related individuals [52]. Similarly, stable 
associations have been reported from insects [50] and 
crustaceans [49] and have also been observed in oyster 
species in the Mediterranean where associations were 
stable even after invasion from the Red Sea [18]. Such 
stable associations harbour an environmental component 
depending on food [49] but also genetic components as 
suggested by similar communities found within mother- 
twin triplets [53]. The fact that the similarity in microbial 
communities correlated with the genetic relatedness of 
the Pacific oyster demonstrated here, further suggests that 
bacterial communities are not only unique to individuals 
but can also assemble according to host genotypes. In 
combination with the lack of significant differentiation of 
community structure between oyster beds this suggests 
that larger scale environmental differences between beds 
may play a limited role when compared to host genotype. 
Furthermore, correlations between genetic microbial com- 
munity distances depended to a large degree on OTUs 
only occurring rarely in the communities (Figure 6). This 



suggests that while abundant taxa may lead a generalist 
life style and are found in the majority of host genotypes, 
rare specialists within the community assemble according 
to host genotypes. An alternative explanation for the for- 
mation of genotype specific microbiome associations is 
vertical inheritance [54,55]. While we cannot rule out this 
possibility for Pacific oysters, the transient nature of the 
genotype specific associations suggests that previously en- 
countered disturbance events should also have led to the 
loss of the inherited genotype-specific microbiota. A re- 
covery of genotype specific associations prior to our ex- 
periment therefore rather suggests an uptake from the 
environment. 

Genetic and phenotypic variation among populations 
in marine invertebrates with type III survival curves can 
occur because of uneven reproductive success among in- 
dividual breeders (i.e., sweepstakes reproductive success 
[56]). Bias in reproductive success between spawning 
events potentially led to the genetic differentiation of the 
investigated oyster beds (Figure 1). Given that we did 
not find patterns of genetic differentiation compatible 
with a stepping stone model or with distance-dependent 
gene flow among oyster beds, a successive formation of 
oyster beds from genetically differentiated spatfall events 
in time is more likely. Successive waves of settlement 
from genetically different broodstocks will also lead to 
structure within beds and increase the genetic diversity 
within populations. 

Sweepstakes reproductive success can also lead to link- 
age disequilibrium, because a reduced effective popula- 
tion size will amplify the effect of genetic drift and thus 
create an overrepresentation of certain allelic combina- 
tions within haplotypes [57]. This also applies to linkage 
disequilibrium between selectively neutral markers (in 
this case microsatellites) and genes with functional rele- 
vance, thus representing potential targets of selection. The 
genetic differentiation that we found between populations 
as well as between individuals should therefore be inter- 
preted as a marker for different spatfall events where vari- 
ation in functional genes (e.g. immune genes) involved in 
microbial colonisation can influence the observed associ- 
ation of host genetics - microbiota relationships. 

Disturbance of microbial communities in oyster gills 

With our parallel tag-sequencing approach we were able 
to describe the microbial communities associated with 
Pacific oyster gill tissue in unprecedented detail, yet the 
38,029 reads used in this analysis were not sufficient to 
capture the total taxonomic richness present in single 
oysters. This represents the typical picture found in mar- 
ine microbial communities in general [20] as well as in 
sediment and open water communities from the same 
habitat [58]. The strongly skewed negative binomial dis- 
tribution of OTUs suggests however that the taxonomic 
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resolution was sufficient to reliably estimate bacterial 
alpha diversity expressed as Shannons H' (Figure 2A). 
Additionally, the parallel characterization of microbial 
diversity in a high number of individuals from different 
oyster beds offers a high level of detail and biological 
replication. Previous studies on the characterisation of 
microbial communities associated with oyster species 
have either been focused on a cultivatable subgroup of 
bacteria [5,59] or used techniques of lower taxonomic 
resolution [18] or coverage [15,16,60] and only very re- 
cently pyrosequencing approaches have been used to 
characterize microbiota of oysters [17]. The gill micro- 
bial communities in our study were dominated by OTUs 
affiliated to the a-proteobacterial genus of Sphingomonas 
sp. The a-proteobacteria are dominant in the open water 
of the Wadden Sea, but rather belong to the SAR11 
group [58]. While we also discovered SARll/Pelagibac- 
ter in ambient samples of DB, its overall abundance was 
very low, indicating that the host associated microbiota 
is not an exact mirror image of the environment. This is 
also reflected in gill associated microbial communities of 
other oyster species that differ more strongly from the 
surrounding sea water than for example gut communi- 
ties [18]. The numerical abundance of a-proteobacteria 
in open water could however partly been attributed to 
PCR bias by preferential amplification of sequences from 
this taxonomic group [61]. The dominant genus de- 
tected, was Sphingomonas which contains opportunistic 
species [62] and can also commonly be found in gill tis- 
sue of European plaice Pleuronectes platessa from the 
same region [38]. It was also abundant on freshly pre- 
pared cod in Iceland [63], indicating that this genus can 
reach high numbers on living hosts but is quickly out- 
competed after the hosts death. 

Dominance of a few closely related OTUs has been re- 
ported for other species of oysters. Zurel et al. [18] for 
example found that between 59 - 79% of OTUs in 
Chama spp. oysters in the Red Sea and the Mediterranean 
belonged to OTUs from the class Oceanospirialles closely 
related to the genera Spongiobacter or Endozoicomonas 
(Hahellaceae), which is known for symbiotic associations. 
While we also observed 47 OTUs from the Oceanospir- 
ialles, these were relatively rare (99 reads in total) and only 
a single OTU was affiliated to the family Hahellaceae. 
Similarly, we only found very few OTUs classified as 
Arcobacter spp. (13 OTUs, 16 reads), which represent a 
major and common component of Chilean oysters Tios- 
trea chilensis [60]. This suggests that oyster microbiomes 
can have similar structures in terms of abundances but 
dominant taxa differ strongly between species, habitats 
and sampled tissues. Under certain environmental condi- 
tions gut communities of other Crassostrea species were 
found to be dominated by Mycoplasma [17], which also 
became dominant in some oysters after disturbance in our 



experiments (Figure 5A). The natural dominance of Myco- 
plasma in oysters from much warmer habitats [17] may 
thus suggest that Mycoplasma represents a temperature 
sensitive part of oyster microbiota and may proliferate 
preferentially at higher temperatures. 

Host stress and abiotic disturbance both could have 
contributed to the major shift in microbial community 
structure (Figure 3). The direction and magnitude of the 
shift was dependent on the initial community compos- 
ition, and although no significant differences were ob- 
served between oyster beds in ambient conditions there 
was some indication for oyster bed specific shifts 
(Figures 3 and 4). The strongest shifts occurred in the 
beds with initially high microbial diversity (OW and 
PK), manifested in a sharp decrease in microbial di- 
versity. In the oyster bed with low diversity on the 
other hand we observed no significant change in bac- 
terial diversity (Figure 2). Our disturbance treatment 
represented a drastic change in abiotic conditions and 
probably resulted in multifaceted host stress re- 
sponses as well as community shifts. We cannot dis- 
entangle what component of stress (food, transfer, or 
heat stress) or microbial community response caused 
the observed shifts. Our aim was however to compare 
the undisturbed natural community to a disturbed 
community in stressed hosts under conditions that 
can facilitate disease outbreaks (i.e., heat waves, food 
depletion, accumulation of waste products). We could 
not observe an overall net increase of obvious patho- 
gen candidates like Vibrio [5,59]. Only OTUs affili- 
ated to Mycoplasma, which can cause disease in 
shellfish [3], showed a strong increase in disturbed 
communities (Figure 4). Mycoplasma were also found 
to dominate microbialcommunities in the gut of Eastern 
oysters Crassostrea virginica [17]. However, since genus 
affiliation will not be sufficient to reliably identify patho- 
genic strains, controlled infection experiments are needed 
to evaluate the true pathogenic potential of the strains 
detected here. Furthermore, since we could neither in- 
voke disease nor observe an increase in the abundance 
or occurrence it seems unlikely that disease agents are 
a constitutive part of the oyster microbiome, suggest- 
ing that disease outbreaks arise from environmental 
sources. 

Mycoplasma was also the taxon that showed the stron- 
gest shift towards a specialist lifestyle (highly abundant 
in few hosts, [46,47], Figure 5 A) and mainly drove the 
trend for higher abundances of specialist taxa in oysters 
exposed to disturbance. This shift towards higher de- 
grees of specialisation also resulted in a positive relation- 
ship between the number of oysters hosting a specific 
OTU (i.e., occupancy) and the mean relative abundance 
of the respective OTU, which was absent from the ambi- 
ent communities (Figure 5A). Such a positive relationship 
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between abundance and occupancy is the null-expectation 
[45] and its absence under ambient conditions can prob- 
ably be attributed to the frequent occurrence of rare taxa 
assembling in a genotype specific manner. On the other 
hand, only a small subset of OTUs shared between treat- 
ments were actually spreading and increasing in abundance 
(mainly Actinobacteria, Sphingomonas and Mycoplasma) 
while others got selectively lost in stressed oysters (mainly 
Flavobacterid). 

Conclusion 

In winter months the microbiome in gill tissue of the in- 
vasive Pacific oyster, Crassostrea gigas, is dominated by 
few highly abundant taxa but show a high taxonomic di- 
versity with many rare taxa supporting previous observa- 
tions from microbial communities in marine sediments 
[20,58]. The (3-diversity of natural, ambient communities 
correlated with individual host relatedness rather than 
with genetic differentiation between oyster beds suggest- 
ing that communities are stable within individuals [18,51] 
and that rare species are associated with genetic differenti- 
ation of the host. This association was lost when the host 
was stressed by our disturbance treatment (Figure 6). The 
subsequent shift in microbial community structure was 
driven by a loss of rare bacterial strains and an increase in 
the abundances of specialist strains. This suggests that 
genotype-specific associations are the result of the overall 
community diversity including rare phylotypes. If it is true 
that the disturbance of ambient host genotype - microbial 
community associations are an important component in 
the defence against infections, it will be very difficult to 
control disease by for example administering probiotics. 
Therefore, monitoring microbial communities during an 
actual infection will be an important future avenue of re- 
search to address the role of genotype specific microbial 
communities for host fitness and to improve our ability to 
predict mass mortality events in benthic populations. 

Availability of supporting data 

Data are available at http://dxdoi.org/10.1594/PANGAEA8198% 
Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

KMW planned the research, performed molecular labwork, and led the 
writing of the manuscript, NV conducted the experimental field and lab 
work, data analyses was done by KMW, HP and AE. All authors read and 
approved the final manuscript. 

Acknowledgements 

We would like to thank three anonymous reviewers for their helpful 
comments and Kevin Schiele for drawing the map in Figure 1. This study 
was financially supported by the Emmy-Noether grant WE4641 -1 of the 
German Research Foundation DFG given to KMW. 



Author details 

^elmholtz Centre for Polar and Marine Research, Alfred Wegener Institute, 
Coastal Ecology, Wadden Sea Station Sylt, Hafenstrasse 43, 25992, List/Sylt, 
Germany, department of Biological Sciences, University of South Carolina, 
715 Sumter Street, Columbia, SC 29205, USA. 3 Benthic Ecology Laboratory, 
IFREMER, DYNECO, BP70, Plouzane 29280, France, laboratory of Aquatic 
Photobiology and Plankton Ecology, Institute of Ecology, University of 
Innsbruck, Technikerstrasse 25, Innsbruck 6020, Austria, department of 
Ecology and Genetcis, Uppsala University, Limnology, Norbyvagen 18D, 
Uppsala 75236, Sweden. 

Received: 4 July 2013 Accepted: 1 November 2013 
Published: 9 November 2013 

References 

1. Harvell CD, Kim K, Burkholder JM, Colwell RR, Epstein PR, Grimes DJ, 
Hofmann EE, Lipp EK, Osterhaus ADME, Overstreet RM, et al: Emerging 
Marine Diseases-Climate Links and Anthropogenic Factors. Science 1999, 
285(5433):1505-1510. 

2. Li Y, Qin J, Abbott C, Li X, Benkendorff K: Synergistic impacts of heat 
shock and spawning on the physiology and immune health of 
Crassostrea gigas: an explanation for summer mortality in Pacific 
oysters. Am J Physiol Regul Integr Comp Physiol 2007, 293(6):R2353-R2362. 

3. Paillard C, Le Roux F, Borrego J: Bacterial disease in marine bivalves, a 
review of recent studies: Trends and evolution. Aquat Living Resour 2004, 
17(4):477-498. 

4. Friedman CS, Estes RM, Stokes NA, Burge CA, Hargove JS, Barber BJ, Elston 
RA, Burreson EM, Reece KS: Herpes virus in juvenile Pacific oysters 
Crassostrea gigas from Tomales Bay, California, coincides with summer 
mortality episodes. Dis Aquat Organ 2005, 63(1)33-41 . 

5. Gamier M, Labreuche Y, Garcia C, Robert A, Nicolas JL: Evidence for the 
involvement of pathogenic bacteria in summer mortalities of the Pacific 
oyster Crassostrea gigas. Microb Ecol 2007, 53(2):1 87-1 96. 

6. Soletchnik P, Lambert C, Costil K: Summer mortality of Crassostrea gigas 
(Thunberg) in relation to environmental rearing conditions. J Shellfish Res 
2005, 24(1 ):1 97-207. 

7. Fleury E, Huvet A: Microarray Analysis Highlights Immune Response of 
Pacific Oysters as a Determinant of Resistance to Summer Mortality. 
Marine Biotech nol 2012, 14(2):203-217. 

8. Samain JF, Degremont L, Soletchnik P, Haure J, Bedier E, Ropert M, Moal J, 
Huvet A, Bacca H, Van Wormhoudt A, et al: Genetically based resistance to 
summer mortality in the Pacific oyster (Crassostrea gigas) and its 
relationship with physiological, immunological characteristics and 
infection processes. Aquaculture 2007, 268(1 -4):227-243. 

9. Wendling CC, Wegner KM: Relative contribution of reproductive 
investment, thermal stress and Vibrio infection to summer mortality 
phenomena in Pacific oytsers. Aquaculture 2013, 412-413:88-96. 

10. Schulenburg H, Kurtz J, Moret Y, Siva-Jothy MT: Ecological immunology. 
Philos Trans R Soc B Biol Sci 2009, 364(1 5 1 3):3-1 4. 

11. Zilber-Rosenberg I, Rosenberg E: Role of microorganisms in the evolution 
of animals and plants: the hologenome theory of evolution. 

FEMS Microbiol Rev 2008, 32(5):723-735. 

12. Dubilier N, Bergin C, Lott C: Symbiotic diversity in marine animals: the art 
of harnessing chemosynthesis. Nat Rev Microbiol 2008, 6(1 0):725-740. 

13. Castro D, Pujalte MJ, Lopez-Cortes L, Garay E, Borrego JJ: Vibrios isolated 
from the cultured manila clam {Ruditapes philippinarum): numerical tax- 
onomy and antibacterial activities. J Appl Microbiol 2002, 93(3):438-447. 

14. Prado S, Romalde JL, Barja JL: Review of probiotics for use in bivalve 
hatcheries. Vet Microbiol 2010, 145(3-4):187-197. 

15. Green TJ, Barnes AC: Bacterial diversity of the digestive gland of Sydney 
rock oysters, Saccostrea glomerata infected with the paramyxean 
parasite, Marteilia sydneyi. J Appl Microbiol 2010, 109(2):61 3-622. 

16. Hernandez-Zarate G, Olmos-Soto J: Identification of bacterial diversity in 
the oyster Crassostrea gigas by fluorescent in situ hybridization and poly- 
merase chain reaction. J Appl Microbiol 2006, 1 00(4):664-672. 

17. King GM, Judd C, Kuske CR, Smith C: Analysis of Stomach and Gut 
Microbiomes of the Eastern Oyster [Crassostrea virginica) from Coastal 
Louisiana, USA. PLoS ONE 2012, 7:12. 

18. Zurel D, Benayahu Y, Or A, Kovacs A, Gophna U: Composition and 
dynamics of the gill microbiota of an invasive Indo-Pacific oyster in the 
eastern Mediterranean Sea. Environ Microbiol 201 1, 13(6):1 467-1 476. 



Wegner et al. BMC Microbiology 2013, 13:252 
http://www.biomedcentral.com/1471-2180/13/252 



Page 12 of 12 



19. Fernandez-Piquer J, Bowman JP, Ross T, Tamplin ML: Molecular analysis of 
the bacterial communities in the live Pacific oyster (Crassostrea gigas) 
and the influence of postharvest temperature on its structure. J Appl 
Microbiol 2012, 1 12(6):1 134-1 143. 

20. Sogin ML, Morrison HG, Huber JA, Mark Welch D, Huse SM, Neal PR, Arrieta 
JM, Herndl GJ: Microbial diversity in the deep sea and the underexplored 
"rare biosphere". Proc Natl Acad Sci U S A 2006, 103(32):1 21 15-12120. 

21. Reise K: Pacific oysters invade mussel beds in the European Wadden Sea. 
Senckenbergiana maritima 1998, 28:167-175. 

22. Buttger H, Nehls G, Witte S: High mortality of Pacific oysters in a cold 
winter in the North-Frisian Wadden Sea. Helgoland Mar Res 201 1, 
65(4):525-532. 

23. Moehler J, Wegner KM, Reise K, Jacobsen S: Invasion genetics of Pacific 
oyster Crassostrea gigas shaped by aquaculture stocking practices. J Sea 
/?es2011 ; 66(3):256-262. 

24. Watermann BT, Herlyn M, Daehne B, Bergmann S, Meemken M, Kolodzey H: 
Pathology and mass mortality of Pacific oysters, Crassostrea gigas 
(Thunberg), in 2005 at the East Frisian coast, Germany. J Fish Dis 2008, 
31:621-630. 

25. Linnenbrink M, Wang J, Hardouin EA, Kuenzel S, Metzler D, Baines JF: 
The role of biogeography in shaping diversity of the intestinal microbiota 
in house mice. Mol Ecol 2013, 22(7):1 904-1 91 6. 

26. Hedgecock D, Pudovkin Al: Sweepstakes reproductive success in highly 
fecund marine fish and shellfish: a review and commentary. Bull Mar Sci 
201 1,87(4):971 -1002. 

27. Oberbeckmann S, Fuchs BM, Meiners M, Wichels A, Wiltshire KH, Gerdts G: 
Seasonal Dynamics and Modeling of a Vibrio Community in Coastal 
Waters of the North Sea. Microb Ecol 2012, 63(3):543-551. 

28. Malham SK, Cotter E, O'Keeffe S, Lynch S, Culloty SC, King JW, Latchford JW, 
Beaumont AR: Summer mortality of the Pacific oyster, Crassostrea gigas, 
in the Irish Sea: The influence of temperature and nutrients on health 
and survival. Aquaculture 2009, 287(1 -2):1 28-1 38. 

29. Li G, Hubert S, Bucklin K, Ribes V, Hedgecock D: Characterization of 79 
microsatellite DNA markers in the Pacific oyster Crassostrea gigas. 
Mol Ecol Notes 2003, 3(2):228-232. 

30. Hubert S, Hedgecock D: Linkage maps of microsatellite DNA markers for 
the pacific oyster Crassostrea gigas. Genetics 2004, 168(1)351-362. 

31. Schuelke M: An economic method for the fluorescent labeling of PCR 
fragments. Nat Biotechnol 2000, 18(2):233-234. 

32. Weir BS, Cockerham CC: Estimating F-Statistics for the Analysis of 
Population- Structure. Evolution 1984, 38(6):1 358-1 370. 

33. Hedgecock D, Li G, Hubert S, Bucklin K, Ribes V: Widespread null alleles 
and poor cross-species amplification of microsatellite DNA loci cloned from 
the Pacific oyster, Crassostrea gigas. J Shellfish Res 2004, 23(2):379-385. 

34. Chapuis MP, Estoup A: Microsatellite null alleles and estimation of 
population differentiation. Mol Biol Evol 2007, 24(3):62 1-631. 

35. Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred 
from metric distances among DNA haplotypes: application to human 
mitochondrial DNA restriction data. Genetics 1992, 131:479-491. 

36. Huber J, Morrison H, Huse S, Neal P, Sogin M, Mark Welch D: Effect of PCR 
amplicon size on assessments of clone library microbial diversity and 
community structure. Environ Microbiol 2009, 11(5):1 292-1 302. 

37. Binladen J, Gilbert MTP, Bollback JP, Panitz F, Bendixen C, Nielsen R, 
Willerslev E: The Use of Coded PCR Primers Enables High-Throughput Se- 
quencing of Multiple Homolog Amplification Products by 454 Parallel 
Sequencing. PLoS ONE 2007, 2(2):e197. 

38. Wegner KM, Shama LNS, Kellnreitner F, Pockberger M: Diversity of immune 
genes and associated gill microbes of European plaice Pleuronectes 
platessa. Estuar Coast Shelf Sci 201 2, 108:87-96. 

39. Quince C, Lanzen A, Davenport RJ, Turnbaugh PJ: Removing Noise From 
Pyrosequenced Amplicons. Bmc Bioinformatics 201 1, 12:38. 

40. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello 
EK, Fierer N, Pena AG, Goodrich JK, Gordon Jl, et al: QIIME allows analysis 
of high-throughput community sequencing data. Nat Methods 2010, 
7(5):335-336. 

41 . Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid 
assignment of rRNA sequences into the new bacterial taxonomy. 

Appl Environ Microbiol 2007, 73(1 6):5261 -5267. 

42. Kunin V, Engelbrektson A, Ochman H, Hugenholtz P: Wrinkles in the rare 
biosphere: pyrosequencing errors can lead to artificial inflation of 
diversity estimates. Environ Microbiol 2010, 12(1):1 18-123. 



43. Austin B, Austin DA: Bacterial Fish Pathogens - Disease of Farmed and Wild 
Fish. 4th edition. Berlin: Springer; 2007. 

44. R Development Core Team: R: A language and environment for statistical 
computing. Vienna: R Foundation for Statistical Computing; 2012. 

45. Gaston KJ, Blackburn TM, Greenwood JJD, Gregory RD, Quinn RM, Lawton 
JH: Abundance-occupancy relationships. J Appl Ecol 2000, 37:39-59. 

46. Barberan A, Bates ST, Casamayor EO, Fierer N: Using network analysis to 
explore co-occurrence patterns in soil microbial communities. ISME J 
2012, 6(2)343-351. 

47. van der Gast CJ, Walker AW, Stressmann FA, Rogers GB, Scott P, Daniels TW, 
Carroll MP, Parkhill J, Bruce KD: Partitioning core and satellite taxa from 
within cystic fibrosis lung bacterial communities. ISME J 201 1 , 5(5):780-791 . 

48. Durban A, Abellan JJ, Jimenez-Hernandez N, Latorre A, Moya A: Daily 
follow-up of bacterial communities in the human gut reveals stable 
composition and host-specific patterns of interaction. FEMS Microbiol Ecol 
2012, 81(2)427-437. 

49. Freese HM, Schink B: Composition and Stability of the Microbial 
Community inside the Digestive Tract of the Aquatic Crustacean 
Daphnia magna. Microb Ecol 201 1, 62(4):882-894. 

50. Robinson CJ, Schloss P, Ramos Y, Raffa K, Handelsman J: Robustness of the 
Bacterial Community in the Cabbage White Butterfly Larval Midgut. 
Microb Ecol 201 0, 59(2):1 99-21 1 . 

51. Vanhoutte T, Huys G, De Brandt E, Swings J: Temporal stability analysis of 
the microbiota in human feces by denaturing gradient gel 
electrophoresis using universal and group-specific 16S rRNA gene 
primers. FEMS Microbiol Ecol 2004, 48(3)437-446. 

52. Lozupone CA, Stombaugh Jl, Gordon Jl, Jansson JK, Knight R: Diversity, 
stability and resilience of the human gut microbiota. Nature 2012, 
489:220-230. 

53. Reyes A, Haynes M, Hanson N, Angly FE, Heath AC, Rohwer F, Gordon Jl: 
Viruses in the faecal microbiota of monozygotic twins and their 
mothers. Nature 2010, 466(7304):334-U381. 

54. Geib SM, Jimenez-Gasco MM, Carlson JE, Tien M, Jabbour R, Hoover K: 
Microbial Community Profiling to Investigate Transmission of Bacteria 
Between Life Stages of the Wood-Boring Beetle, Anoplophora glabripennis. 
Microb Ecol 2009, 58(1 ):1 99-21 1 . 

55. Inoue R, Ushida K: Vertical and horizontal transmission of intestinal 
commensal bacteria in the rat model. FEMS Microbiol Ecol 2003, 
46(2):213-219. 

56. Li G, Hedgecock D: Genetic heterogeneity, detected by PCR-SSCP, among 
samples of larval Pacific oysters (Crassostrea gigas) supports the hypoth- 
esis of large variance in reproductive success. Can J Fish Aquat Sci 1998, 
55(4):1 025-1 033. 

57. Nei M, Li W-H: Linkage disequilibrium in subdivided populations. Genetics 
1973,75:213-219. 

58. Gobet A, Boer SI, Huse SM, van Beusekom JEE, Quince C, Sogin ML, Boetius 
A, Ramette A: Diversity and dynamics of rare and of resident bacterial 
populations in coastal sands. ISME J 2012, 6(3)542-553. 

59. Lacoste A, Jalabert F, Malham S, Cueff A, Gelebart F, Cordevant C, Lange M, 
Poulet SA: A Vibrio splendidus strain is associated with summer mortality 
of juvenile oysters Crassostrea gigas in the Bay of Morlaix (North 
Brittany, France). Dis Aquat Organ 2001, 46:139-145. 

60. Romero J, Garcia-Varela M, Laclette JP, Espejo RT: Bacterial 16S rRNA gene 
analysis revealed that bacteria related to Arcobacter spp. constitute an 
abundant and common component of the oyster microbiota (Tiostrea 
chilensis). Microb Ecol 2002, 44(4)365-371. 

61 . Gonzalez JM, Moran MA: Numerical dominance of a group of marine 
bacteria in the alpha-subclass of the class Proteobacteria in coastal sea- 
water. Appl Environ Microbiol 1997, 63(9361410)4237-4242. 

62. Piccini C, Conde D, Alonso C, Sommaruga R, Pernthaler J: Blooms of single 
bacterial species in a coastal lagoon of the southwestern Atlantic Ocean. 
Appl Environ Microbiol 2006, 72(10):6560-6568. 

63. Reynisson E, Lauzon HL, Magnusson H, Jonsdottir R, Olafsdotir G, 
Marteinsson V, Hreggvidsson GO: Bacterial composition and succession 
during storage of North-Atlantic cod {Gadus morhua) at superchilled 
temperatures. BMC Microbiol 2009, 9(1 9961 579):250. 



doi:1 0.1 1 86/1 471 -21 80-1 3-252 

Cite this article as: Wegner et al.: Disturbance induced decoupling 
between host genetics and composition of the associated microbiome. 

BMC Microbiology 201 3 1 3:252. 



